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Abstract 

A comprehensive model for the photon number fluctuations and the current noise in quantum cascade 
lasers is presented. It is shown that the photon intensity noise in quantum cascade lasers exhibits little 
amplitude squeezing even when noise in the drive current is suppressed below the shot noise value. This 
is in contrast to interband semiconductor diode lasers in which the laser intensity noise can be squeezed 
well below the shot noise limit by high impedance suppression of fluctuations in the drive current. The 
theoretical model presented in this paper self-consistently accounts for the suppression of current noise 
in electron transport in multiple quantum well structures due to various electronic correlations. The 
nature of these electronic correlations is discussed. Mechanisms responsible for the reduced photon num- 
ber squeezing in intersubband lasers are elucidated. Scaling of the laser intensity noise and the current 
noise with the number of cascaded gain stages is also described. Direct current modulation response of 
quantum cascade lasers is also studied, and it is shown that contrary to the predictions in the literature 
of terahertz modulation bandwidth for these lasers, bandwidth of almost all quantum cascade lasers that 
have been reported in the literature is limited by the inverse photon lifetime inside the laser cavity to 
tens of gigahertz. 



78.67.De, 72.70.+m, 73.21.Fg 



Contents 



1 Introduction 



2 Outline 



3 Rate Equations and Steady State Solutions 



3.1 Steady State Solutions 



3.1.1 Below Threshold 



3.1.2 Above Threshold 



4 

7 
7 



4 Theoretical Model for Noise and Fluctuations 



4.1 Linearized Langevin Rate Equations for Electron and Photon densities 10 

4.2 Linearized Electron Transport, Coulomb Correlations and Noise 10 

4.3 Displacement Currents 14 

4.4 Differential Resistance 15 



4.5 Electron Transport in the Supcrlatticc Injector 16 

4.6 Biasing Electrical Circuit 16 



5 Solution of the Coupled Equations 



18 

5.1 Current Modulation Response) 18 



5.2 Laser Intensity Noise and Current Noise 21 



6 Current Noise: Results and Discussion 



23 

3.1 Circuit Models for the Current Noisd 23 



6.2 Spectral Density and Fano Factor of the Current Noise 23 



6.2.1 Effect of Small Differential Impedance of a Single Gain Stage 27 

6.2.2 Effect of Electronic Correlations 27 

30 

30 



6.3 Scaling of the Current Noise with the Number of Cascade Stages 

6.4 Spectral density of the Current Noise in the External Circuit . . 



7 


Photon Noise: Results and Discussion 






7.1 Spectral Density and Fano Factor of the Laser Intensity Noise 






7.2 Scaling of the Laser Intensity Noise with the Number of Cascade Stages 




7.3 Effect of Multiple Longitudinal Modes on the Measured Intensity Noise 





8 Conclusion 



A Correlations Among the Langevin Noise Sources 



B Differential Resistance of a QCL 



C Elements of Matrix D 



30 

30 
36 
36 

37 

37 

38 

38 

38 

40 

40 

41 

G.l Modulation Response 42 

G.2 Differential Resistance 43 

G.3 Differential Impedance 44 

G.4 Current Noise] 44 

G.5 Suppression of the Current Noise by Large External Impedance] 45 

G.6 Intensity Noise 45 



D Important Elements of Matrix D 



E Noise Spectral Densities and Fano Factors 



F Current Noise and the Equivalent Circuit Model; 



G Noise Model for Semiconductor Quantum Well Diode Lasers 



List of Figures 



Superlattice quantum cascade laser 



Multiple quantum well quantum cascade laser 



Energy subbands of the three levels of the gain stage. Most favored electronic transitions by 



optical phonon emission are almost horizontal in the E(k) — k space 



\i Electron densities in level 3 and level 2 of the gain stage, and the output power per facet are 



plotted as a function of the bias current. There is a discontinuity in the rate of increase of the 



electron densities with the bias current at threshold 



5 Charge densities associated with the electron densities Snjnj, Sn^, Sn2, and 5n\ are shown 



The electron charge densities are imaged on the positively charged ionized dopants present in 
the superlattice injector of the subsequent stage 



|6 Differential Resistance of a QCL is shown as a function of the bias current. The experimen- 



tally measured discontinuity in the differential resistance at threshold is about 0.3 ft. The 



theoretical model reproduces the discontinuity exactly. The experimental data is taken from 



7 Circuits used for biasing QCLs 



Ref. [11 



S Thevenin equivalent circuit model indicating the distinction between 5I ext (uj) and SI meas (uj) 



9 Direct current modulation response (20 log 10 [|_ff (wlj/LHYO)]]) of QCLs is plotted as a function 



of the frequency for different bias currents. Modulation response shown in the figure has been 



normalized w.r.t its value at zero frequency. For the device parameters see Table 1 



10 Absolute value of the QCL impedance Z(uj) is plotted as a function of the frequency for differ- 



ent bias currents. The peaks in the values of Z(uj) are not because of relaxation oscillations 



since the modulation response of the QCL is over-damped, but because the smallest zero of 



Z(lu) is smaller than its smallest pole 



11 Circuit model for the current fluctuations 



12 A simplified circuit model for the current fluctuations 



13 Spectral density Ki(uo) of the current noise is plotted as a function of the frequency. The 



noise spectral density has been normalized w.r.t. its value at zero frequency 



14 Low frequency spectral density Ki(lu = 0) of the current noise is plotted as a function of the 



bias current 



15 Fano Factor for the low frequency current fluctuations is plotted as a function of the bias 



current 



16 N times the Fano Factor for the low frequency current fluctuations is plotted as a function ol 



the bias current. 



17 Fano Factor for the (low frequency) noise in the laser intensity is plotted as a function of the 



bias current 



18 Low frequency Relative Intensity Noise (RIN) is plotted as a function of the bias current 



Very small amount of squeezing (less than 0.4 dB) is exhibited at high bias levels even when 



the circuit current fluctuations arc suppressed with a 50 Q impedance 



19 Low frequency Spectral density Kp(uj = 0) of the laser intensity fluctuations is plotted as a 



function of the bias current. When Z s is large (50 Q) small amounts of squeezing is possible 



at high bias levels 



20 Relative Intensity Noise (RIN) for the QCL with improved r\ Q is shown. Only around 1 dB ol 



squeezing is possible at high bias levels and when Z s = 50 f2 



21 Relative Intensity Noise (RIN) is plotted as a function of the frequency for different bias 



currents (Z s = O). At high frequencies the RIN reaches the shot noise limit 



22 Fano Factor for the laser intensity noise is plotted as a function of the frequency (Z s = £1) 



The intensity noise at frequencies higher than the inverse of the photon lifetime in the cavity 



is dominated by the photon partition noise at the output facet 



23 Ratio of the RIN obtained by ignoring the term containing the current fluctuations in Equation 



(96|) to the actual RIN is plotted as a function of the bias current for different values of the 
impedance Z s . The current fluctuations are suppressed when Z s is large and the error incurred 
in calculating RIN is, therefore, small. 



24 Active region of a semiconductor quantum well diode laser. 



25 Circuit model for the current fluctutions in semiconductor diode lasers 



1 Introduction 



Unipolar quantum cascade lasers (QCLs) utilizing intersubband transitions to generate photons have become 
important sources of light in the mid- infrared region (5 /^m - 15 /xm). In this paper a comprehensive model 
for the photon intensity fluctuations in QCLs is reported for the first time. Current noise associated with 
electron transport through the active regions is also studied and its effect on the photon intensity noise is 
evaluated. 

QCLs are different from interband semiconductor diode lasers in three important ways which can have a 
significant impact on their noise properties: 

1. Electron transport in QCLs takes place by tunneling between states in adjacent quantum wells. It is 
well known that electronic correlations in resonant tunneling in quantum well structures can suppress 
(or enhance) current noise by providing a negative (or positive) feedback (l], |^, Sj. High impedance 
suppression of the current noise in semiconductor diode lasers results in light output with squeezed 
photon number fluctuations . It is therefore intriguing whether suppression of the current noise can 
also lead to squeezing in QCLs. Any model for the photon noise in QCLs must take into account these 
electronic correlations self-consistcntly. 

2. In diode lasers the carrier density in the energy level involved in the lasing action does not increase 
beyond its threshold value and, therefore, the noise contributed by the non-radiative recombination 
and generation processes also remains unchanged beyond threshold. In QCLs the electron densities in 
the upper and lower lasing states do not clamp at threshold, and keep increasing when the bias current 
is increased beyond threshold. As a result, non-radiative processes contribute significantly to photon 
noise even at high bias currents. 

3. Since all the gain sections in a QCL are connected electrically and optically, electron density fluctuations 
and photon emission events in different gain sections become correlated. The effect of these correlations 
on the photon noise in interband cascade lasers has already been discussed in detail || |(| , and it is the 
aim of this paper to investigate the role of these correlations in QCLs. 

2 Outline 

In section || the non-linear rate equations for the electron and photon densities in QCLs are presented. The 
steady state solution of these rate equations below and above threshold are described. In section || the 
non-linear rate equations are linearized to obtain Langevin rate equations for the fluctuations in the electron 
and photon densities. Electron transport in the multiple quantum well structure of QCLs is discussed in 
detail, and a self-consistent model for the fluctuations in the electron charge densities and the electron 
current density is presented. It is shown that a self-consistent description of the fluctuations in the charge 
and current densities can be carried out in terms of a few device parameters. Langevin noise sources are 
also used to model the noise associated with electron transport by tunneling. Section |^ is the main part 
of this paper. In section || the set of coupled linearized Langevin rate equations for the fluctuations in 
the electron densities in different levels of all the cascaded gain stages and the fluctuations in the photon 
density are solved under the constraints imposed by the biasing electrical circuit. In addition, the direct 
current modulation response of QCLs is also evaluated and the maximum possible modulation bandwidth 
is discussed. The analytical and numerical results on the current noise and the photon noise in QCLs are 
presented and discussed in section |^ and section ^, respectively. In these sections the results obtained are 
compared with the current and photon noise in interband semiconductor diode lasers. Readers not familiar 
with the results on the current and photon noise in diode lasers are encouraged to read Appendix |G| in which 
a detailed model for the noise in diode lasers is presented. 

3 Rate Equations and Steady State Solutions 

Many different types of QCL structures have been reported in the literature [Q, |[ ^, [h], [ll], [l2| [b| [ll| 
|l6| , [l?], |l^, [HJ. Almost all of these QCL structures can be classified into two categories: 



1. Superlattice QCLs in which the gain stage consists of a superlattice structure and the photons are emit- 
ted when the electrons make transitions between two minibands of this superlattice. These minibands 
are actually clusters of closely spaced energy levels (Fig. [jl6| [l7| [l8[ |l9) . 

2. Multiple quantum well QCLs in which the gain stage consists of multiple quantum wells (typically two 
or three) and the radiative electronic transitions occur between two discrete energy levels (Fig. |2|) jjj 

In both types of QCLs, two successive gain stages are separated by a superlattice structure known as the 
injector. The superlattice injector has a mini-gap which prevents the electrons from tunneling out into the 
injector from the upper energy level(s) of the previous gain stage and, therefore, increases the radiative 
efficiency. Electrons from the lower energy level(s) of a gain stage can tunnel into the injector, and the 
injector injects these electrons into the upper energy level(s) of the next gain stage. 

In this paper, photon noise and current noise in only multiple quantum well QCLs is discussed. However, 
the methods described here can easily be used to study noise in QCLs with superlattice gain stages. A single 
stage of a multiple quantum well QCL structure is shown in Fig. || jll) . The operation of the QCL can be 
described as follows. Electrons tunnel from the energy states in the superlattice injector into level 3 of the 
gain stage. Photons are emitted when electrons make radiative transitions from level 3 to level 2. Transitions 
from level 2 to level 1 occur primarily by emission of optical phonons. Electrons leave the gain stage from 
level 1 by tunneling out into the injector of the next stage. In addition, electrons also make non-radiative 
transitions from level 3 to levels 2 and 1. The non-linear rate equations for the electron and photon densities 
can be written as, 

dr 'k- J L- R 32 (4,ni) ' RsM,4) - Vv g 9(4,4) (Sp + (1) 



dt q 

= r 32 (4,4) - R 2 i(4M) + r% sK 4) (s P + (2) 



^ = r 31 (4, 4) + r 2 i(4,4) - — (3) 

dt q 

^ = £r S9( „i„j) (s , + ^)-| ,4, 

Pout = Vo hv WLS p (5) 

In the above equations, n J k is the electron density (cm -2 ) in the fcth energy level of the jth gain stage. J( n 
and J 3 out are the electron current densities (cm -2 ) tunneling into level 3 and tunneling out of level 1 of the 
jth gain stage, respectively. Only in steady state Jf n equals J 3 out - S p is the photon density (cm -2 ) inside 
the optical cavity. S p is equal to the total number of photons inside the cavity divided by the width W 
and the length L of the cavity. v g is the group velocity of the lasing mode and g is the optical gain (cm -1 ) 
contributed by a single gain stage. is the mode confinement factor for the jth gain stage. N is the total 
number of cascaded gain stages. R32 is the net transition rate from level 3 to level 2 through non-radiative 
processes and spontaneous emission into the non-lasing modes. Similarly, R31 and R21 are the net transition 
rates from level 3 and level 2 into level 1, respectively. n sp is the spontaneous emission factor (20). P ou t is 
the output power from the laser. r/ is the power output coupling efficiency and t p is the photon lifetime 
inside the cavity. The expression for t p is, 



1 

— = Vg(cti + a m ) = v g 



1 1 ( 1 



(G) 



where cti is the internal loss of the cavity, a m is the loss from the cavity facets, and r\ and r 2 are the facet 
reflectivities. The power output coupling efficiency n from the facet with reflectivity r\ is, 
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Figure 1: Superlattice quantum cascade laser. 
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Figure 2: Multiple quantum well quantum cascade laser. 



For simplicity it is assumed that all the gain stages have the same mode confinement factor, i.e. r j = T 
for all j. This assumption is valid if all the cascaded gain stages are located close to the peak of the transverse 
profile of the optical mode where the field strength varies slowly. Even for QCLs with large number of gain 
stages numerical simulations show that corrections to the solution obtained by assuming all r j to be equal 
are small. Under this assumption, the steady state electron densities nL are the same in all the gain stages, 
and the index j may be suppressed. 



3.1 Steady State Solutions 
3.1.1 Below Threshold 

The steady state solution to the rate equations can be found by setting all the time derivatives equal to zero, 
and putting Ji n — J. Below threshold, steady state carrier densities can be found by putting S p — 0, and 
solving the equations (the index j has been suppressed), 

R32(n 3 ,n 2 ) + R 3 i(n 3 ,ni) = — (8) 

fl32(n 3 ,n 2 ) = -R2i(n 2 ,ni) (9) 
The third equation can be obtained by realizing that J out is also a function of n\ , 

Jout{ni) = J (10) 

To proceed further, analytical expressions for the transition rates are required. These transition rates can 
be approximated as, 



R32(n 3 ,n 2 ) = — 


(11) 


T~32 


R3i(n 3 ,ni) = — 


(12) 


T31 


-K2i( n 2,m) = — 


(13) 


T21 


Jout{n\) ^ ni 


(14) 


q T out 



The rationale for the approximations in Equations (|ll|)-([l^) is that optical phonons are largely responsible 
for intersubband transitions. As shown in Fig. [| optical phonon mediated intersubband transitions that 
are almost horizontal in E(k)-k plane are more likely to occur Therefore, the transitions rates from an 
upper to a lower subband are not much affected by the electron density in the lower subband, as long as the 
electron density in the lower subband is small. More complicated expressions for these transition rates, such 
as, 

R q k{n q , n k ) = - — (15) 

may be used if necessary. 

The expression for J out in Equation (^) does not depend upon on the electron density in the injector since 
electrons in the injector states are assumed to relax very quickly into the ground state of the injector which 



is spatially localized near the next gain stage. Using Equations (|TT|)-([l4j) in Equations (g)-(10), expressions 
for the carrier densities can be obtained as a function of the current density, 

J T32T 3 1 

n 3 = (16) 

q T32 + T31 

n 2 = ■ (17) 

q 7-32 + T31 

m = —T out (18) 
q 




3.1.2 Above Threshold 

Above threshold, the gain is clamped to a value determined by equating the gain with the loss, 

N 1 

^2 r % ff(«3' n i) = NTv g ff("3, n 2 ) = — (19) 

3 = 1 1 

For perfectly parabolic subbands, the expression for the gain may be approximated as, 

9{n 3 ,n 2 ) = a(n 3 - n 2 ) (20) 

where a is the differential gain. The carrier and photon densities above threshold can be obtained by solving 
the equations, 

#3i(n3,ni) +i?2i(n 2 ,ni) = — + — = - (21) 

T31 T 2 i q 

q Tout q 

NTv g a(n 3 - n 2 ) = — (23) 



which results in, 



n 3 = - T21 T31 + ( 1 ^ T 3i ^ 

gr 2 i+T3i \NTvgaTpJ t 2 i + T31 

n2 = l^U^„( 1 (25) 
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Figure 4: Electron densities in level 3 and level 2 of the gain stage, and the output power per facet are 
plotted as a function of the bias current. There is a discontinuity in the rate of increase of the electron 
densities with the bias current at threshold. 



Pout=VoVr—N(I-I th ) (28) 

q 

where the threshold current density Jth and the radiative efficiency r\ r are, 

Jth = ~xrp~ (— + —) 7j ( 29 ) 

NTVgdTp \T 32 T 31 J (1-T21/T32J 
Vr = (l - ^) ^ 

The radiative efficiency rj r for a QCL is defined as that fraction of the total number of electrons injected 
into each gain stage which contribute to photon emission. Equation ( p8|) shows that the slope efficiency of 
the laser scales linearly with the number of gain stages N. 



Equations (|24[) and (|25|) show that above threshold, even though the gain is clamped to its threshold 
value, the electron densities keep increasing with the bias current. This is in contrast to semiconductor diode 
lasers in which the carrier densities do not increase beyond their threshold values. As a result, an increase 
in the injected current density in QCLs does not only lead to an increase in the photon emission rate but 
it also leads to an increase in the rate of non-radiative transitions. Therefore, QCLs tend to have radiative 
efficiencies i] r significantly smaller than unity. Fig. ^ shows the electron densities n 3 and n 2 plotted as a 
function of the bias current. The values of the various device parameters used in generating Fig. |] belong 
to the QCL reported in pH , and these values are given in Table |]. Fig. |5] shows that the rate of change of 
electron densities in levels 3 and 2 with the bias current exhibits discontinuities at threshold. This can be 
confirmed by comparing Equations ( p4| ) and ( ps| ) with Equations ( |l6| ) and (|f7|). As will be shown later in 
this paper, these discontinuities in the rate of increase of electron densities with the bias current result in a 
discontinuity in the value of the differential resistance of the laser at threshold. 



4 Theoretical Model for Noise and Fluctuations 



The model for the noise presented in this paper consists of a set of coupled self-consistent Langevin rate 
equations for the fluctuations in the electron density in different energy levels of a gain stage. Fluctuations 



in the electron density are caused by radiative and non-radiative scattering processes, electron tunneling 
processes and also by fluctuations in the current injected into the gain stage. Fluctuations in the current 
are a relaxational response to electron scattering and tunneling events occurring inside all the gain stages 
of the QCL, and they are also caused by sources external to the laser which include thermal noise sources 
associated with circuit resistances. Photon density fluctuations are also modeled by Langevin rate equations. 
Electron density fluctuations in different gain stages are all coupled to the photon density fluctuations and 
also to the fluctuations in the current which flows through all the gain stages connected in series. The system 
of equations obtained this way can easily be solved analytically or numerically to give the spectral density 
of the photon number fluctuations and the current fluctuations. The methods described in this paper can 
be used to study a variety of QCLs that have been reported in the literature. 



4.1 Linearized Langevin Rate Equations for Electron and Photon densities 

The non-linear rate equations can be linearized around any bias point to obtain rate equations for the 
fluctuations. Linearized Langevin rate equations for these fluctuations are, 
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Equations (|ll|)-([T3|) have been used above for approximating the transition rates R q k- /32, /31, and /21 are 
Langevin sources which model the noise associated with the non-radiative intersubband transitions and also 
the radiative transitions into the non-lasing modes. Jrn and fas are Langevin sources which model the 
noise in photon emission and absorption from the lasing mode. Fl and F a describe the noise associated 
with photon loss from the cavity ]20| . All the Langevin noise sources have a white spectral density and their 
correlations can be found by the methods described in Ref. ^0|. All the non-zero correlations among the 
noise sources are given in Appendix [A]. 



4.2 Linearized Electron Transport, Coulomb Correlations and Noise 

In order to determine S J{ n and J 3 out the electron transport through the active region needs to be looked at in 
detail. Self-consistent modeling of electron transport in multiple quantum well structures poses a significant 
challenge, and the steady state current-voltage characteristics of QCLs are difficult to compute accurately. 
In this paper a slightly different approach has been adopted which is more useful for the problem under 
consideration. A self-consistent model for the fluctuations in the electron current density and the electron 
charge density is presented. It is shown that self-consistent analysis of current density and charge density 
fluctuations can be carried out in terms of only a few device parameters. The values of these parameters can 
cither be determined experimentally or computed theoretically from more detailed self-consistent transport 
models. The method used in this paper to estimate the value of each parameter will be discussed when we 
compare the theoretical model with the experimental results. 



The expression for the direct sequential tunneling current density from the injector state into level 3 of 
the gain stage can be written as p2j, 
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where the forward and backward components of the injection current are, 
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p is the coupling constant, and its related to the transmission probability. Ei n j(k) and E 3 (k) are the 
energies of electrons in the injector state and level 3 of the gain stage, respectively. A(E) is a normalized 
lorentzian function with FWHM equal to the broadening of the energy levels, and f(E — is the Fermi- 
Dirac distribution function with a chemical potential /i. Expressions similar to Equation ( p6|) can also be 
written for the phonon assisted tunneling current density. The analysis presented in this paper is independent 
of the specific nature of the electron tunneling mechanisms. In what follows, and E 3 will stand for 



Einj(k = 0) and E 3 (k = 0), respectively. The tunneling current in Equation (36) depends upon the following 
three quantities, 



• The difference {^i n j 
state. 



Einj) between the injector chemical potential and the energy of the injector 



• The difference (fi 3 — E 3 ) between the chemical potential and the energy of level 3 of the gain stage. 

• The relative difference (Ei n j — E 3 ) between the energies of the injector state and level 3 of the gain 
stage. 

The current can change if the number of electrons in the injector level or in level 3 of the gain stage changes. 
The current can also change if the energy of the injector level shifts with respect to the energy of level 3. 
5Ji n can be written as, 
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fin is a Langevin noise source which models the noise in electron tunneling. Noise in electron transport 
by sequential tunneling in multiple quantum well structures can be described with Langevin noise sources. 
In Refs. ^ m the current noise in double barrier resonant tunneling structures is evaluated using classical 
discrete master equations. Under suitable conditions a discrete master equation may be converted into a 
Fokker-Planck equation, and if the fluctuations are relatively small a Fokker-Planck equation can be linearized 
around a stable steady state solution (see Ref. p3| for details) . Langevin rate equations can be used in place 




Figure 5: Charge densities associated with the electron densities 5m n j, 6713, 5ri2, and 5ni are shown. The 
electron charge densities are imaged on the positively charged ionized dopants present in the superlattice 
injector of the subsequent stage. 



of linearized Fokker-Planck equations since the two formalisms are equivalent. [ p3| . It can be shown that 
Langevin rate equations yield results identical to those presented in Refs. B [| for the current noise in double 
barrier resonant tunneling devices p4}| . A linearized analysis based on Langevin rate equations may become 
invalid for highly non-linear devices. The correlation function for the noise source fi n is, 

WL{fi n (t) ff n {t')) = ~ (Jin- forward + Jin-backward) SjgS(t -t') (41) 

* — XinS jq S(t-t') (42) 

The factor Xin relates the sum of the forward and backward tunneling currents to their difference which is the 
total injection current Ji n . At low temperatures Xin is expected to be close to unity since Pauli's exclusion 
would restrict the available phase space for the backward tunneling current ||. For the same reason Xin is 
expected to be close to unity for large values of the injection current Ji n . At high temperatures and small 
values of the injection current Xin can be larger than unity. 

Although Equation ( f40| ) for the change in current density is derived for the case of direct sequential 
tunneling, it also holds for the case of phonon assisted tunneling. Even if the energy distribution of electrons 
inside each energy level in steady state were not a Fermi-Dirac distribution, Equation ( p0| ) would still hold. 

It is assumed that the superlattice injector is doped in regions not close to the gain stage. Electric held 
lines from electron density fluctuations Sn J in p 5n J 3 , Sir^, and 5n\ are imaged on the ionized dopants in the 
injector layer of the (j + l)th gain stage, as shown in Fig. ||. Therefore, the fluctuation SV^ in the potential 
difference across the jth gain stage can be written as, 

5 y i = ^<L + ^A + ^A + ^ (43) 
(^inj ^3 ^2 C/l 

Cinj, C3, C2, and C\ are capacitances which determine the incremental change in the potential difference 
across a gain stage with changes in the electron densities in different energy levels. Using first order quantum 
mechanical perturbation theory, SE J in ^ — 6E3 can be related to the fluctuation in the average potential 



difference between the injector ievei and level 3 of the gain stage. The fluctuation in the average potential 
difference between these two levels can also be expressed in terms of capacitances. The expression for 
SEj n j — SE 3 therefore becomes, 
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In Equation ( |4^ ) Gi n , t%. 



, and £3 are given by, 
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(47) 



More generally, there may be more than just one energy level in the injector from which electrons get 
injected into level 3 of the gain stage. Equation (45) can be modified appropriately to take into account the 
contributions from all the energy levels inside the injector. However, if the values of tj n are roughly the same 
for all such states in the injector then the final form of Equation (|45| ) will remain unchanged, but Srii n j will 
then represent the total electron density in all the injector states. 

Similarly, the fluctuation S J 3 out in the tunneling current density from level 1 of the gain stage into the 
injector is given by the expression, 
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In a well designed QCL the backward tunneling current from the injector of the next stage into level 1 of 
the gain stage is small, and Xout is expected to be close to unity. E' 3 n - is the energy of the injector level of 
the next stage into which electrons tunnel from level 1 of the gain stage. SE^ — SE 3 n p as before, can be 
expressed in terms of capacitances, 



5E j -SE' j - g2(5 "'"J 

orj l 0Cj inj — r"' 



q 2 Sn 3 3 



q 2 8n 3 2 



q 2 Sn 3 1 

~cf~ 



(52) 



Using Equation (p3|) and Equation rt5 
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where tout and G ou t are, 

1 1 o J ou t 5J out . . 

7 — = — x ' G °«* = IjT^ (54) 

tout q oni 6(Ei - E' inj ) 

In Equation ( ji8| ) it has been assumed that electrons in the injector relax into the ground state state of the 
injector sufficiently fast so that electron occupation in the injector levels do not effect the electron escape 
rate out of level 1 of the gain stage. 

Note that Gi n and G ou t can be positive or negative depending upon the relative alignment of the energy 
levels Ei n j and E3 in the steady state. The scheme used in deriving Equation (fl6|) and Equation (|53| ) is 
fairly general and can be used to derive self-consistent linearized transport equations for a variety of QCL 
structures. Approximations can be made to simplify Equation ([46| ) and Equation ([33]). Expression for SJ? n 
can also be written as, 
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For the sake of economy of notation new parameters have been introduced in the above equation, 

1 1 j Gi n 1 1 j Gi„ 1 _ Gi„ 1 _ Gi„ 

Tin t-in G^j T3 t3 G3 T2 Gj 7i G{ 

Simple electrostatic arguments can be used to show that T2 and n will be large, and can be assumed to be 
infinite. 

The injector is assumed to have a large number of closely spaced energy levels. J ou t is, therefore, largely 
insensitive to the relative shifts in E\ and E( •. This implies that terms containing G ou t in the expression 
for 5J 3 out may be neglected, and the simplified expression for 5J out becomes, 

^ = —S4 + fLt (58) 
q Tout 

where T ou t is just t ou t- Equations (^6|) and ( |58| ) show that in addition to the parameters given in the electron 
and photon density rate equations (Equations (|3ll)-(j34|)), the paramters necessary for describing electron 
transport through the gain stage are Ci n j , G3 , C2 , C\ , Ti„ , T3, T2, and t\ . 

4.3 Displacement Currents 

The noise current SJ ex t, which flows in the external circuit, is not equal to <5Jf n or 8J 3 out . 5J ex t also includes 
displacement currents and is given by the expression, 

d Sti 3 

SJe Xt = 6Jl l + q ^f^ (59) 

Since all the gain stages are connected electrically in series, the same current 5J ex t flows through all the 
gain stages. The second term on the right hand side of Equation (|59| ) is the contribution to SJ ext from 



displacement currents. Differentiating both sides of Equation (43) with respect to time and rearranging 
yields, 

dnL l _ dSV^_^ C mj dSnj 
q dt ~ inj dt ^ q C k dt 

Using the above equation the expression for SJ ex t becomes, 
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Equation (pTI) follows from Equation (pOj) by using the particle number conservation equation, 
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Equation Q60J) and Equation (|61j) satisfy the Ramo-Shockley theorem ]25| 

4.4 Differential Resistance 

Below threshold, the total differential resistance -fi^ of all the gain stages can be calculated by substituting 
Equations 



|§) in Equation (|56|), 
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Above threshold, the differential resistance can be computed by using Equations 



(62) 
(63) 

i) with Equation 



N 



WLC, 



N 



in] 

J_C ri 

Tin 



1 



WLC, 



Co 

■ (1 



1 Ci„j 
Tin C3 
_J_\ _T 2 lT 3 l 
T2 



T21 



T31 
l) 



T21 + r 3 i 

1 Ci n j 
r in Cl 



(64) 
(65) 



">3 



Expressions for the parameters 9$, 9' 3 , 9 2 , 9' 2 , and 9\ are given in Appendix [B]. Notice the similarity between 
Equations (|6^) and ( |65| ) , and Equation ( |195| ) for the differential resistance of interband semiconductor diode 
lasers given in Appendix |G]. Unlike the active regions of diode lasers, the active regions of unipolar QCLs 
are not charge neutral, and as a result various capacitances appear in the expression for the differential 
resistance of QCLs. 

The discontinuity ARd in the differential resistance at threshold for an TV stage QCL is, 
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The Incremental change in the potential drop across a gain stage is related to the incremental changes 
in electron densities through Equation (ff^). Therefore, the discontinuity in the differential resistance at 
threshold results from the discontinuity at threshold in the rate of change of electron densities in levels 
3 and 2 of the gain stage with the bias current. Fig. || shows the calculated and measured differential 
resistance of a QCL as a function of the bias current. The experimental data is taken from Ref. |pH . The 
values of the various device parameters are given in Table |. Values of T21, t 3 i and T32 are taken from 
Rcf. pj . Values of all the capacitances given in Table [j] are estimated from the structure of the QCL 
described in Ref. pi]] . Values of Tj n , T3, and T out are estimated from Equations ( f47| ) and (54). The total 
resistance of the ohmic contacts and the superlattice injectors is assumed to be approximately 0.3 at 
threshold. The experimentally observed discontinuity in the differential resistance at threshold is exactly 
reproduced in the calculated results without the use of any fitting parameters. This agreement suggests that 
the self-consistent model for the linearized electron transport presented in this paper adequately captures 
the essential ingredients. 



10r 



^ 9- 

a 

_i 

o 

2 n 



v 
o 
c 

*" 5 ■ 

.52 

CO 

o 

a: 
re 

c 

a 

a 



Q 1 



experiment 
theory 



discontinuity (AR d ) in the differential 
resistance at threshold 




0.2 0.4 0.6 0.3 

Bias Current (A) 



Figure 6: Differential Resistance of a QCL is shown as a function of the bias current. The experimentally 
measured discontinuity in the differential resistance at threshold is about 0.3 fl. The theoretical model 
reproduces the discontinuity exactly. The experimental data is taken from Ref. PH . 



Diode lasers also exhibit a discontinuity in the differential resistance at threshold. As shown in Ap- 
pendix |G|, the discontinuity in the differential resistance of diode lasers at threshold is K B T/qI t h times a 
factor of the order of unity, which can be compared with the more complicated expression given in Equation 
© for QCLs. 

4.5 Electron Transport in the Superlattice Injector 

In this paper no attention has been given to modeling the electron transport through the superlattice 
injector. In the absence of any bias current, the energy levels in the injector are not suitably aligned 
to facilitate electron transport, and the resistance of the injector region is large. As the bias current is 
gradually increased electrons pile up in different quantum wells until their presence modifies the potential 
profile and aligns the energy levels such that the electron current can flow. Once the injector has been turned 
on in this fashion, the differential resistance of the injector region is negligible, and the only bottleneck for 
electron transport is the gain stage. As a result of the small differential resistance of the injector region, 
any current noise originating in the injector region will not couple well into the external circuit. Therefore, 
electron transport in the injector region may be ignored when modeling noise. If necessary, the impedance 
of the superlattice injectors can be modeled with a lumped element, and the current noise generated inside 
the injector regions can be modeled with a voltage source in series (or a current source in parallel) with that 
element, as shown below. A detailed discussion of the current noise in superlattice structures is beyond the 
scope of this paper. 

4.6 Biasing Electrical Circuit 

Two electrical circuits for biasing QCLs are shown in Fig. [?]. In circuit A the QCL, with an impedance 
Z(ui), is biased with a voltage source V s in series with an impedance Z s (ui). The thermal noise originating 
in the impedance Z s (u>) is modeled by adding a voltage noise source SV S . The differential impedance of the 
superlattice injector and the current noise generated by the injector can also be modeled with an impedance 
and a voltage noise source in series (or a current noise source in parallel) with that impedance. For the 
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Figure 7: Circuits used for biasing QCLs. 



sake of economy of notation it will be assumed that the impedance Z a (u) represents not just an external 
circuit impedance but the Thevenin equivalent impedance of the superlattice injectors, device ohmic contacts, 
external circuit resistances, and device and circuit parasitics, and the voltage noise source 8V S represents the 
Thevenin equivalent of their individual noise sources. Only the gain stages inside the QCL are not included 
within Z s {ui) and they are represented by the impedance Z[uS). However, Z(u>) will be loosely referred to as 
the impedance of the QCL. The current noise generated by the gain stages can also be modeled by adding 
a current noise source in parallel with Z{ui), as shown in later sections. 

Direct current modulation of the QCL can be achieved by adding an RF voltage source in series with V s , 
and this RF voltage source can also be represented by the voltage source 8V S . From the context it will be 
clear whether 8V S represents a RF signal source or a noise source. 

In circuit B the QCL is biased with a current source in series with an ideal inductor, and it is also 
capacitively coupled to a voltage source 8V S ' with a series impedance Z' s {oj). If at frequencies of interest 
the inductor and the coupling capacitor are almost open and short, respectively, then this circuit is also 
equivalent to circuit A. Therefore, in this paper only circuit A will be considered. In circuit A the current 
SI ext can be expressed as, 

N 



8V s {uj) - ^2SV j (u 



8I ext {u) = 8J ext (w) WL = -L± (68) 

Z s (ui) 

It is important to note here that SI ex t(u>) may not be the noise current which would be measured in an 
experiment. For example, suppose that the QCL has a parasitic capacitance C in parallel with the actual 
device, as shown in Fig. ||. The QCL is driven with a series resistor R Q and a noise voltage source 8V (ui) 
representing the thermal noise in the resistor R Q . Fig. |^ shows the distinction between the noise current 
SI ex t(oj) defined in Equation (|68|), and the noise current 8I meas (oS) that would be measured in an experiment. 
Notice that the Thevenin equivalent impedance Z s (uj) is a parallel combination of the resistance R and the 
capacitance C . Z s (u>) and 8V s (ui) are, 
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Figure 8: Thevenin equivalent circuit model indicating the distinction between 5I ex t{u) and <5/ meas (a;). 



and the relation between 8I ex t{u>) and 5I mea s{^>) is, 

5I ext (u) = 



8I r , 



(l+jw£(w)C ) 



(70) 



Choosing to define Z s (lo) this way helps in formulating a noise model that is independent of the specific 
nature of the device parasitics. 



5 Solution of the Coupled Equations 

5.1 Current Modulation Response 

In this section the response 8P ou t(Ld) / 8I ex t(uj) of QCLs to external sinusoidal current modulation 8I ex t(u) 
is determined []. The frequency dependence of the photon noise spectral density of semiconductor lasers is 
directly related to the frequency dependence of the current modulation response. It is therefore instructive 
to look at the modulation response of QCLs. The modulation response can be found by solving Equations 
(|3l|)-(|35|), together with Equations ( |56|) and (|58|) , and setting all the noise sources equal to zero. The external 
circuit constraints expressed in Equations (|60[) and ([38]) must also be enforced. Equations (|3l|)-(|33|) for each 
gain stage are coupled to the same set of equations for all the other gain stages through Equations ( |34| ) 
and (]60|). Such a large system of coupled equations can be solved only numerically. A numerical approach, 
although simple to implement, is not very instructive. With the approximation that all gain stages have the 
same confinement factor T, a significant portion of the work can be done analytically. This approach will be 
followed in this paper. All equations, unless stated otherwise, will be expressed in the frequency domain. 

The relationship between the current density 8J ex t(oj), which flows in the external circuit, and the total 
potential drop 5V(to) across all the gain section can be obtained by using Equation ( |56| ) in Equation (|60|), 
and summing over the index j, 
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: It is assumed that P a ut(t) = Pout + Rea,l{SP out {ui) e^*} and I ext (t) = I + Keal{&I ext {u) e^*} 



The following new symbols have been introduced in Equation (|7l|), 
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Using Equations (p6|), (|58| ) and (|7l| ) in Equations (|3l|)-(|34[), summing over the index j, and arranging the 
resulting equations in a matrix form gives, 
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The coefficients of the matrix D can be found from Equations (|3l|)-(|34|), and they are given in Appendix |c]. 
The solution of Equation d72) can be written as, 
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The coefficients of the matrix D 1 are given in Appendix [D[ Equation (73) can be used in Equation (|7l]) to 
calculate the total impedance Z(uj) of all the gain stages, 
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Z(cj = 0) is just the differential resistance Rd of the QCL given earlier in Equations ( |62|) and (|64|). Finally, 
from Equation (35) and Equation (|73j), the current modulation response can be written as, 
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Above threshold, an analytical approximation for the modulation response can be found in the limit WTj„ <C 1 
and {t2 , n} — * 00 (Appendix |d|), 
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(7"21 + T31) 

Expression for the current modulation response function H(u>) is given in Appendix [D| In Equation (76) r st 
is the differential lifetime associated with stimulated and spontaneous photon emission into the lasing mode, 
and is given by the relation, 

For frequencies less than the inverse of the smallest time constant of the QCL (which is usually either T2i 
or T s t), Equation (|7^) can be further simplified, and put in the standard form used for semiconductor diode 
lasers (see Appendix and Ref. p0|), 
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where r\ r is the radiative efficiency defined in Equation (p0|), and the relaxation oscillation frequency u>r and 
the damping constant 7 are, 
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The damping constant 7 can be related to lur, 
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The JsT-factor describes the damping of the QCL modulation response at high photon densities. 7 G has a 
weak dependence on the photon density through r st , and it approaches 1/t p at large photon densities. 

If the condition ujr < j/V2 is satisfied then Equation ( |79| ) describes a second order over-damped system. 
For many QCLs that have been reported in the literature this condition holds true above threshold. Using the 
values of device parameters from Table |], ujr and 7 can be calculated. If we assume that the output power 
of the laser is around 150 mW, then from Equations and ((78j) t p and r s t are approximately 7 ps and 2.8 
ps, respectively. The resulting value of 7 is more than three times larger than that of lor. The internal time 
constants in QCLs are usually smaller than the photon lifetime t p , and therefore the modulation response 
of QCLs is over-damped. In contrast, the current modulation response of semiconductor diode lasers is 
under-damped, and becomes over-damped only at very large bias currents when r s t becomes small po| . 

For QCLs the 3 dB frequency, which is defined to be the frequency at which the square modulus of the 
laser modulation response becomes one half of its value at zero frequency, can be found from the simplified 
expression for the modulation response in Equation (|79|), 
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As the photon density inside the laser cavity increases the 3 dB frequency also increases but it asymptotically 
approaches an upper limit W3dB|max- This maximum attainable 3 dB bandwidth can be calculated from 
Equation (|S5|), and it comes out to be, 

w 3dB I max ~ (86) 

Tp 

As long as the computed value of W3dB | max is less than 1/Yj n , I/V21, and l/r s t, the approximations made in 
deriving Equation (|7^) are justified. Equation (^) confirms the intuition that a laser cannot be modulated 
much faster than the inverse of the photon lifetime inside the laser cavity. As shown in Appendix ^[ in 
diode lasers the value of ui 3 dB | max equals \/2/t p . The difference of a factor of \/2 comes from the fact that 
in diode lasers the modulation response is under-damped (see Appendix |c|) . 

As in diode lasers, the photon lifetime imposes a fundamental limit on how fast QCLs can be modulated. 
It is not uncommon to find predictions of THz modulation bandwidths for QCLs in literature [p6| . However, 
for all the QCLs reported in the literature so far, the photon lifetime is the longest of all the time constants 
and it is the dominant factor that would limit the modulation bandwidth of these QCLs to tens of GHz 
instead of THz. 




Figure 9: Direct current modulation response (20 log 10 [|ff (u;)|/|if (0)|]) of QCLs is plotted as a function of 
the frequency for different bias currents. Modulation response shown in the figure has been normalized w.r.t 
its value at zero frequency. For the device parameters see Table [P. 



Fig. |9| shows the calculated modulation response of a QCL as a function of the frequency for different 
values of the bias current. The values of the different parameters of the QCL are taken from Ref. Ill] , 
and are given in Table In the numerical calculations values of all the device time constants (except r sf ) 
were assumed to be independent of the bias. Fig. ^ shows that at low bias currents the 3 dB frequency 
increases with the bias current, and at high bias currents the 3 dB frequency saturates to a value which is 
well approximated by 1/(2ttt p ) = 21 GHz. The analysis carried out in this paper does not take into account 
device heating which may also be important in limiting the modulation bandwidth of QCLs at large biases. 

Fig. [To] shows the impedance Z(u) of the QCL plotted as a function of the frequency for different bias 
currents. The peaks in the values of Z(oj) are not due to relaxation oscillations since, as already pointed 
out earlier, the modulation response of the QCL is over-damped. The peaks are due to the fact that the 
smallest zero of Z(ui) is at a frequency which is smaller than the frequency of its smallest pole. Impedance 
measurements can therefore provide valuable information about the time scales associated with electron 
dynamics in QCLs. 

5.2 Laser Intensity Noise and Current Noise 

In this section the current noise in the external circuit and the intensity noise in the output power from 
QCLs is calculated. In the Langevin equation formalism noise is added through the Langevin noise sources 
which were introduced in Equations (|3~f|)-(|35|) and also in Equations ( |56| ) and (|5^). In addition, any noise 
originating in the external circuit and in the superlattice injectors can also contribute to the current noise 
and the photon noise, and as already explained earlier, this noise can be represented by the voltage source 
SV S . In this paper it is assumed that SV S represents the thermal noise originating in the series impedance 
Z s (uj), and its correlation function is, 

(5V s (uj) 5V s {u>')) = 2K B T Resd{Z s (uj)} 2n 5(uj - uj') (87) 

By assuming the above correlation function for the noise source SV S we are ignoring any noise which may be 
contributed by the superlattice injectors. 
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Figure 10: Absolute value of the QCL impedance Z(ui) is plotted as a function of the frequency for different 
bias currents. The peaks in the values of Z(u>) are not because of relaxation oscillations, since the modulation 
response of the QCL is over-damped, but because the smallest zero of Z(uj) is smaller than its smallest pole. 



Including the Langevin noise sources Equation (72) can be written as, 
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The expressions for the noise sources Fx, F2, F3, and F4 are, 
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The solution of Equation (p 



F 4 (w) 
can be written as, 



JV 

-E 

3=1 



F L (m) 
N 



' 5Ni(u) ' 




- d^h ■ 


5N 2 (uj) 












SSp(uj) 




. D 43 H 



N 8J ext {uo) 



Ef=i ^V) 
Ef=i E>3» 

Eli DiV) 



(89) 

(90) 
(91) 
(92) 



(93) 



8J ex t(u) in Equation ( J93| ) is an unknown. Using Equation (|5E|) in Equation (|60|), summing over the index j, 
and making use of Equation yields, 
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where Fj„ = S^Li //«• Substituting the value of 6V(ui) from Equation ( |68| ) in Equation (|94|), we get the 
final expression for the current fluctuations SI ex t(io) in the external circuit, 
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The fluctuation SP out (u>) in the output power can be obtained by substituting Equation (|9^) in Equation 
(|93|), and using Equation (|35|), 
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6 Current Noise: Results and Discussion 

6.1 Circuit Models for the Current Noise 

Equation ( |95|) for the current fluctuations in the external circuit can be written as, 
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Expression for SP(u>) is given in Appendix ^. Observing that Z(lu)/N is the differential impedance of a 
single gain stage, a circuit model for the current fluctuations can be constructed by attaching a current noise 
source 5P(u>) in parallel with the jth gain stage, as shown in Fig. [fl]. Current noise sources belonging to 
two different gain stages are not independent but are positively correlated, as shown in Appendix [f| This is 
because carrier density fluctuations in all the gain stages interact with the same optical field. 

A simplified circuit model for the current noise, more relevant from the point of view of experiments, 
is shown in Fig. |l2|, where a single current noise source 51{us) has been added in parallel with all the gain 
stages of the QCL. Expression for 6I(u>) is, 
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Equation ( p7| ) shows that the current noise SI (u>) is equal to the current noise 5I ex t(io) in the external circuit 
if 5V s (lo) and Z s (co) are both zero. This is also obvious from Fig. [l^. The characteristics of the noise source 
8I{lo) are explored next. 

6.2 Spectral Density and Fano Factor of the Current Noise 

The spectral density Ki(lu) of the noise source 6 1 (to) can be calculated from Equation (|95|). Most of the 
numerical results presented in this paper, unless stated otherwise, are for the QCL described in Ref. [ pd| . 
The device parameters for this QCL are given in the Table [l| In the numerical calculations values of all the 
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Figure 11: Circuit model for the current fluctuations. 
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Figure 12: A simplified circuit model for the current fluctuations. 
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Figure 13: Spectral density Ki(u>) of the current noise is plotted as a function of the frequency. The noise 
spectral density has been normalized w.r.t. its value at zero frequency. 



device time constants (except r st ) were assumed to be independent of bias. The values of Xin an d Xout were 
assumed to be unity (see the discussion in section 4.2). Fig. |l^ shows the frequency dependence of Kj(uj) 
for different values of the bias current. As expected, Kj{u) rolls over near the 3 dB frequency (w3dB) for the 
laser modulation response. Fig. [15] shows the Fano Factor (Appendix |E|) for the low frequency fluctuations 
of the current noise source SI(oj) as a function of the bias current. Near the laser threshold the fluctuations 
in the electron densities inside the gain stages become large, and consequently, the current noise is also large. 
Away from the laser threshold the current noise is suppressed far below the shot noise value. 

For frequencies less than wadB, analytical expression for Ki(u>) can be found using the expressions for 
the elements of the matrix D _1 given in Appendix 
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Expressions for the parameters #3, 9' 3 , 9 2 , 9' 2 and 9\ are given in Appendix [B], The expression for Kj{uS) 
above threshold is valid provided, 
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It is insightful to compare the expression for the current noise in Equation (|99) 

interband semiconductor diode lasers. Using the model presented in Appendix |G one gets for diode lasers 
(see Appendix ^ for details), 
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Figure 15: Fano Factor for the low frequency current fluctuations is plotted as a function of the bias current. 



where r/i is the current injection efficiency, 9' is a number of the order of unity, and 9 is much less than 
unity (see Appendix |G|) . Comparing Equations (99) and (101) one can see that below threshold and also 



much above threshold (when r st — > and / 3> Ith) the current noise approaches the shot noise value in 
diode lasers, whereas in QCLs the current noise can be suppressed much below the shot noise value. The 
mechanisms responsible for the suppression of the current noise in QCLs are discussed below. 

6.2.1 Effect of Small Differential Impedance of a Single Gain Stage 

The total differential impedance of all the gain stages in a iV-stage QCL is larger than the differential 
impedance of a single gain stage by a factor of N. This reduces the total noise power of the current 
fluctuations by a factor of N, and therefore Ki(u>) has an explicit l/N dependence in Equation (|99|). 

6.2.2 Effect of Electronic Correlations 

The expression for the current fluctuations 81 \ui) given in Equation (^8j), for frequencies less than W3dB 7 can 
also be written as, 



NSI{uj) _ ^ 
qWL ~ 



(102) 



Equation (102), which is almost identical to Equation (202) given in Appendix |g] for semiconductor diode 



lasers, shows that fluctuations in the electron density in different levels of the gain stage causes fluctuations 
in the current. The sign of the current fluctuations is such as to restore the electron density to its average 
value, thus providing a negative feedback. The physical mechanisms responsible for this negative feedback 
are discussed below. On one hand these electronic correlations suppress the current noise associated with 
electron injection into the gain stage by providing negative feedback, and on the other hand they are also 
responsible for generating current noise in response to electron density fluctuations caused by noise sources 
internal to the gain stage. Various physical mechanisms included in our model which contribute to these 
electronic correlations are described below: 

1. Coulomb Correlations: If the electron density changes in any level of the gain stage then the electrostatic 
potential energy of level 3 also changes because of coulomb interactions. As a result, the energy level 
separation 5Ei n j — 8E3 also changes, and consequently the total electron current from the injector into 
the gain stage also changes. Usually QCLs are not biased in the negative differential regime and the 
value of the conductance Gi„, given by Equation (fl7|), is positive. Therefore, the change in the current 
will be such as to restore the electron density in the levels of the gain stage to its steady state value. 
Coulomb correlations provide negative feedback to regulate electron density fluctuations. If a QCL is 
biased in the negative differential regime, in which the coulomb correlations provide positive feedback 
(negative Gi n ), the fluctuations may increase substantially and the linearized noise analysis presented 
in this paper may not be applicable. In our model the effect of coulomb correlations was introduced 
through the parameters Gi n /C' k in Equation (fl6|). 

2. Pauli's Exclusion and Backward Tunneling Current: If the electron density increases in level 3 of the 
gain stage then this reduces the phase space available for additional electrons to tunnel into level 3 from 
the injector due to Pauli's exclusion, and consequently the forward tunneling current from the injector 
into level 3 decreases from its average value. In addition, an increase in electron density in level 3 also 
increases the backward tunneling current from level 3 into the injector and this also reduces the net 
current from the injector into level 3 (recall from Equation ([36]) that the net current is the difference 
of the forward and backward tunneling currents). In our model both these effects were introduced 
through the parameter t^ 1 in Equation ([hi]). We remind the readers that later in Equations ( |55| ) and 
( pq ) t^ 1 and Gin/C^ were absorbed in the definition of 73 , and Gm/C^ and Gi n /C' 3 were relabeled as 
T2 and t{~ , respectively. Therefore, coulomb correlations, Pauli's exclusion and backward tunneling 



current account for the presence of the terms 8n J k {uj) j Vfc in Equation (102) 



3. Injector Electron Density Response: Here we explain the presence of the terms (Ci n j /ri n CV-) 8n k {oj) 



in Equation (102). Recall that the current fluctuations 8I(uj) can be evaluated by looking at the 
current fluctuations 8I ext (uj) in the external circuit when Z s (u>) is zero, and all external voltage sources 
are incrementally shorted, and the sum of the fluctuations in voltage across all the gain stages (i.e. 



Ej=i SV^) is, therefore, also zero. Under these conditions the relationship between the fluctuations in 
the carrier densities, expressed earlier in Equation ([43|), becomes, 
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Equation (103) can be used to write Equation (|102|) as 
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Equation (104) shows that the current fluctuations are proportional to the total fluctuations in the 
electron density in the injector states of all the stages. Since J2j=i SV^ui) — 0, a net increase in the 
electron density in different levels of all the gain stages must result in a net decrease of the electron 
density in all the injector states, and consequently, the current being injected into the gain stages must 
also decrease. This effect is captured through the terms (Ci n j/Ti n Ck)Sn J k (uj) appearing in Equation 



(102) 



As a result of the electronic correlations described above the current noise associated with electron injec- 
tion into the gain stages, which is represented in our model through the noise sources /iL(w)> is suppressed. 
Electron density fluctuations caused by sources internal to the gain stage contribute more strongly towards 
the current fluctuations because of the same correlations. To see this in a more transparent fashion it is best 
to write Equation ( 102 ) in terms of all the noise sources. Below threshold Equation (102) becomes, 
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Above threshold 5I(u>) is, 
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Note that the strength of the electronic c orre lations dep ends on the values of the parameters #3, 9' 3 , 6*2, 6* 2 
and 0\ (Appendix |b|). From Equations (105) and (106) it is clear that larger values of these parameters 
will result in stronger electronic correlations, larger suppression of the current noise associated with electron 
injection into the gain stage, and also larger contribution to the current noise from the noise sources internal 
to the gain stage. The reader is encouraged to compare Equat ions ( |105| ) and (106) with the corresponding 
expressions for semiconductor diode lasers given in Equations ( ^03) and (§0j) in Appendix 

A quantitative measure of the role played by the electronic correlations in suppressing the current noise 
can be obtained by multiplying the Fano Factor of the current noise by N. It has been mentioned earlier 
that a factor of 1 /N appears in Equation (|9^) because the total differential impedance of all the gain stages 
is larger than the differential impedance of a single gain stage by a factor of N . Therefore, multiplying the 
current noise Fano Factor by N removes this explicit 1/N dependence in the current noise, and the resulting 
expression can only be less than unity because of electronic correlations. Fig. |l6| shows the current noise 
Fano Factor from Fig. [l^ multiplied by N. Below threshold and much above threshold N times the current 
noise Fano Factor is less than 0.5. This implies that electronic correlations are responsible for suppressing 
the current noise by a factor greater than 2. From Equation (pa), expression for the current noise Fano 
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Figure 16: N times the Fano Factor for the low frequency current fluctuations is plotted as a function of 
the bias current. 



Factor Fj can be written as, 
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For semiconductor diode lasers, using Equation (101), one gets, 
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At frequencies much higher than the inverse of the smallest time constant of the QCL the current noise 
5I(uj) is just the capacitive response to the various electronic transitions which occur inside the gain stages. 
In the limit w — > oo, Ki{u>) is given by the expression, 
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Above threshold an extra term, 
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is added to the above equation to account for the stimulated transitions. Semiconductor diode lasers on the 
other hand are charge neutral. Therefore, in the limit ui — > oo, the current noise in diode lasers is just the 
noise associated with carrier injection into the active region (see Appendix p|), 
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6.3 Scaling of the Current Noise with the Number of Cascade Stages 

In QCLs spectral density Kj(lu) of the current noise obeys a simple scaling relation with respect to the 
number of cascaded gain stages N, and this relation can be determined from Equation (pq), 



N 2 K r (u;,I/I th ,N) = N^K^J/I^N') 



(112) 



According to the above equation, the spectral density of the current noise, when expressed as a function of 
J/Jt/i, scales as 1/N 2 . This scaling relation for Kj(uj) holds for all frequencies provided that the transition 
rates Rjk{rij, nk) and the material gain 5(713, 712) are linear functions of the electron densities and the total 
mode confinement factor also scales linearly with the number of cascade stages N. 



6.4 Spectral density of the Current Noise in the External Circuit 

The quantity that can be measured experimentally is the spectral density of the current noise SI ex t(to) which 
flows in the external circuit. Spectral density of the current noise in the external circuit provides a valuable 
tool for studying the high speed dynamics of QCLs. When Z 8 (u>) 7^ SI, which is usually the case, then 
Ki cTt (ui) is not the same as Kj(uj). Expression for if/ cxt (w) follows from Equation (971), 
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Equation ( |113j ) shows that in the presence of a large impedance Z s (uj) the current fluctuations in the external 
circuit are suppressed. The total differential impedance of a QCL is usually less than 1 f2. Therefore, for even 
a moderately large impedance Z s {lo) the current noise in the external circuit is dominated by the thermal 
noise from the impedance Z s (u). Experimental measurement of the current noise would, therefore, require 
a relatively sensitive measurement scheme. High impedance suppression of the current noise in the external 
circuit can influence the laser intensity noise, as shown in Equation (|96| ) and in more detail in the following 
section. 



7 Photon Noise: Results and Discussion 

7.1 Spectral Density and Fano Factor of the Laser Intensity Noise 

The spectral density Kp(lu) of the laser intensity noise can be calculated from Equation (|96|). The Fano 
Factor for the low frequency fluctuations in the laser output power is plotted as a function of the bias current 
in Fig. [It]. The Relative Intensity Noise (RIN) is plotted in Fig. [ll| In each figure the respective shot noise 
limit is also shown. It is assumed that the light coming out from both the facets of the laser is collected 
before the noise is evaluated. This is equivalent to assuming that the output coupling efficiency i] , defined 
earlier in Equation (Q), is, 

Vo = °2 (114) 
(a m + cti ) 

In practice this can be achieved by HR/AR coating the laser facets so that that most of the light comes 
out from only one facet of the laser. When the value of the external impedance Z s is 17, the photon noise 
remains above the shot noise limit. Even at high bias levels, no amplitude squeezing is observed desp ite the 
fact that the current noise is suppressed much below the shot noise value as shown earlier in Fig. [l^. When 




Figure 17: Fano Factor for the (low frequency) noise in the laser intensity is plotted as a function of the 
bias current. 




Figure 18: Low frequency Relative Intensity Noise (RIN) is plotted as a function of the bias current. Very 
small amount of squeezing (less than 0.4 dB) is exhibited at high bias levels even when the circuit current 
fluctuations are suppressed with a 50 O impedance. 
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Figure 19: Low frequency Spectral density Kp(uj = 0) of the laser intensity fluctuations is plotted as a 
function of the bias current. When Z s is large (50 f2) small amounts of squeezing is possible at high bias 
levels. 



Z s = 50 f2, and the current noise in the circuit is further suppressed, a very small amount of squeezing is 
observed at high bias levels (less than 0.2 dB at I = 101th)- This is in sharp contrast to the characteristics 
of interband semiconductor diode lasers which exhibit substantial amount of amplitude squeezing (around 
3-5 dB) when the noise in the circuit current is suppressed far below the shot noise value p7[ |. The absence 
of any significant amount of squeezing in QCLs can be attributed to the following factors, 

1. Small Output Coupling Efficiency: Most of the QCLs reported in literature have relatively long cavity 
lengths and large values for internal waveguide losses. The output coupling efficiency rj a is thus small. 
The small output coupling efficiency significantly reduces the amount of squeezing that can be achieved. 
This, however, does not impose a fundamental limit on the amount of squeezing achievable in QCLs, 
but only a practical limit. 

2. Non-Radiative Transitions: As shown earlier, in QCLs above threshold the electron densities in different 
energy levels of a gain stage do not remain fixed at their threshold values. The electron densities keep 
increasing when the bias current is increased beyond threshold. As a result, the contribution of non- 
radiative electronic transitions to the photon noise also keeps increasing with the bias current. This is 
very different from what happens in semiconductor diode lasers in which above threshold the carrier 
densities in the quantum wells, and the non-radiative recombination rates, remain clamped to their 
threshold values. Amplitude squeezing is therefore expected to be less in QCLs than in semiconductor 
diode lasers, even if the output coupling efficiency of QCLs is improved. The QCL, whose characteristics 
are shown in Fig. ^ and Fig. [S], has a 3 mm long cavity, a waveguide loss of 11 cm -1 , and an output 
coupling efficiency of only 28%. Consider a QCL with a 500 /im long cavity, a waveguide loss of 5 
cm -1 , and an output coupling efficiency of 84%. The values of all the other parameters of this QCL 
are identical to those given in Table [p. Fig. ^ shows the relative intensity noise when this QCL is 
driven with a 50 fl resistor in series. Only about 1 dB of squeezing is observed even at very high bias 
levels (/ ss 101th)- The amount of squeezing achievable in QCLs can be increased by decreasing the 
lifetime t 2 i of electrons in level 2 of the gain stage. This will reduce the rate of increase of the electron 
density in levels 3 and 2 with the bias current above threshold, inrease the radiative efficiency rj r , and 




reduce the contribution of non-radiative electronic transitions to the photon noise. 



Fig. ^l] shows the RIN as a function of the frequency for different values of the bias current assuming 
Z s (lu) = 0. The RIN also rolls over at the frequency 0J3 dB- Fig. shows that the Fano Factor for the laser 
intensity noise as a function of the frequency. As in all other lasers, at frequencies much higher than the 
inverse of the photon lifetime inside the cavity, the RIN is dominated by the noise from photon partition at 
the output facet. Therefore, 
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In this paper careful attention has been given to modeling the current fluctuations in the external circuit. 
The question arises if such detailed modeling of the current fluctuations is necessary for calculating the 
photon intensity noise. In Equation ( p)6| ) the current noise SI e xt{u) is included in the first term on the right 
hand side. It should be noted here that the first and the second term on the right hand side in Equation 
( p6| ) are correlated, and the spectral density of the photon intensity noise cannot be obtained by a simple 
addition of the spectral densities of these two terms. In Fig. ^3| the ratio of the spectral density of the photon 
intensity noise obtained by ignoring the term containing SI ext (u>) in Equation (|9^) to the actual spectral 
density of the photon intensity noise is plotted as a function of the bias current for different values of the 
impedance Z s . When the laser is biased a little above threshold the fluctuations in the current are large, 
and the error incurred by ignoring the term containing 5I ex t(ui) in Equation (Q) is also large. Also, when 
Z s is much larger than the total differential impedance of the QCL, the current fluctuations in the circuit 
are suppressed, and the term containing 5I ex t(u>) can be ignored in Equation (j96|). 

For large Z s (u>), using the expressions for the elements of the matrix D 1 in Appendix [j^, analytical 
expression for the spectral density of the low frequency fluctuations in the laser intensity noise can be 
obtained, 
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Figure 21: Relative Intensity Noise (RIN) is plotted as a function of the frequency for different bias currents 
(Z s = Q). At high frequencies the RIN reaches the shot noise limit. 
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Figure 22: Fano Factor for the laser intensity noise is plotted as a function of the frequency (Z s =0O). 
The intensity noise at frequencies higher than the inverse of the photon lifetime in the cavity is dominated 
by the photon partition noise at the output facet. 
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Figure 23: Ratio of the RIN obtained by ignoring the term containing the current fluctuations in Equation 
( |96| ) to the actual RIN is plotted as a function of the bias current for different values of the impedance 
Z s . The current fluctuations are suppressed when Z s is large and the error incurred in calculating RIN is, 
therefore, small. 
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rj r in the above Equation is the radiative efficiency defined in Equation (pOj). The constants a and (3 are, 
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The above expression for Kp{ui) is valid for frequencies smaller than W3dB, and when the laser is biased 
above threshold when the conditions given by Equation ( |100| ) are satisfied. The expression given in Equation 
(117) is almost identical to the expression for Kp{uS) for semiconductor diode lasers (when the later are also 
biased with a high impedance current source). Using the model presented in Appendix |G] one gets for diode 
lasers (see Appendix ^| for details), 
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where r\i is the current injection efficiency into the quantum wells ]20fl . Comparing Equations (116) and 
(120), it can be seen that the noise associated with the non-radiative transitions appears in a slightly more 
complex form in Equation (116). The noise contribution from the non-radiative transitions in diode lasers 



is constant above threshold and it has been expressed in terms of the threshold current in Equation (120). 
In QCLs, as already pointed out earlier, the noise contribution from the non-radiative transitions keeps 
increasing with bias beyond the laser threshold. Since only a fraction rj r of the electrons injected in level 
3 of the gain stage end up producing photons, a multiplicative factor rfc appears with the noise associated 
with the transitions from level 3 to level 1. A fraction 1 — r\ r of the vacancies left by removing electrons from 
level 2 get filled by radiative transitions from level 3 to level 2, and therefore a factor (1 — i] r ) 2 appears with 
the noise associated with the transitions from level 2 to level 1. All the electrons taken out of level 2 and 
injected into level 3 will end up producing photons since, 

(1 - r, r ) + T) r = 1 (121) 

and, therefore, the noise contribution from the non-radiative transitions from level 3 to level 2 has no 



multiplicative factor in Equation (116). The noise associated with the electron transitions from level 1 
into the injector of the next stage does not directly contribute to photon noise at low frequencies. These 
transitions contribute to the current noise in the external circuit, which can in turn contribute to the photon 



noise. But we have assumed in Equation (116) that Z s (u>) is large and the current fluctuations are suppressed. 
Similarly the noise associated with the electron transitions from the injector into level 3 of the gain stage is 
also suppressed at low frequencies when Z s (oj) is large. 

In diode lasers, since the current injection efficiency rji is less than unity, the partition noise associated 
with carrier leakage from the SCH region contributes a term to the intensity noise which increases linearly 



with the bias current even beyond the laser threshold, as shown in Equation (120). Since r\% is usually close 
to unity in well designed diode lasers ]2C| ], the contribution of this term to the intensity noise is small. 

The Fano factors for the laser intensity noise much above threshold (when r st — ► and I ^> I t h) in QCLs 



and diode lasers can be calculated from Equations (117) and (]120|). 



r , _> l-VoVr + ^o— — ^ ( — ) QCLS , , 

[w < u 3 dB , 1 > 1th) 1 _ Diode Lasers 



Equation (122) gives the maximum photon number squeezing possible in QCLs and in diode lasers. In diode 
lasers both rji and r\ Q can be larger than 0.85 ^o). Therefore, the intensity noise in diode lasers can be 
suppressed more than 5 dB below the shot noise value. For the QCL whose parameters are listed in Table |], 
r/ r and rj have the values 0.66 and 0.28, respectively, and, consequently, the maximum possible squeezing is 
only 0.6 dB. 

7.2 Scaling of the Laser Intensity Noise with the Number of Cascade Stages 

In QCLs K p (uj) obeys a simple scaling relation with respect to the number of cascaded stages N, and this 
relation can easily be deduced from Equation (|96|), 

K P (uj,I/I th ,N)=Kp(uj,I/I th ,N') (123) 

According to Equation ( p.23[ ), the spectral density of the photon noise, when expressed as a function of I/Ith, 
is independent of the value of N. The scaling relation for K p (u) holds for all frequencies and when Z s (uo) is 
very large or when Z s (u>) = O, provided that the transition rates R(rij, n q ) and the material gain g(n 3 , n 2 ) 
are linear functions of electron densities and the total mode confinement factor also scales linearly with the 
number of cascaded stages N. In Ref. |^8) it is shown that the total mode confinement factor scales with 
the number of cascaded stages according to the expression ~ er/(0.019iV), which is almost linear in N for 
N < 40. 

7.3 Effect of Multiple Longitudinal Modes on the Measured Intensity Noise 

Most QCLs reported in literature lase with multiple longitudinal modes. Although the intensity noise of 
each longitudinal mode can be large, the intensity noise of all the modes taken together is expected to be 
adequately described by the single mode analysis carried out in this paper. This is because the intensity 
noise in different lasing modes is negatively correlated, as it is in the case of semiconductor diode lasers p9[ . 
However, this demands that in experiments designed to measure the intensity noise attention must also be 
paid to optimizing the light collection efficiency such that photons are collected from all the lasing modes, 



otherwise intensity noise in excess of that described by Equation (116) can be introduced 



8 Conclusion 



A comprehensive model for treating noise and fluctuations in intersubband quantum cascade lasers has been 
presented. The current noise exhibited by QCLs is much below the shot noise value. Suppression of the 
current noise in QCLs is largely due to the small differential resistance of individual gain stages compared 
to the total differential resistance of all the cascaded gain stages. In addition, electronic correlations also 
tend to suppress the current noise. However, unlike semiconductor diode lasers, current noise suppression 
does not lead to significant photon number squeezing in QCLs. In QCLs the contribution to the photon 
intensity noise coming from the non-radiative electronic transitions keeps increasing with bias beyond the 
laser threshold, and this reduces the amount of photon number squeezing achievable in QCLs compared to 
semiconductor diode lasers. In addition, poor laser output coupling efficiencies and the lack of availability of 
photon detectors with good quantum efficiencies in the mid-infrared wavelength region make the experimental 
observation of photon number squeezing in QCLs difficult. 

The current modulation response of QCLs has also been investigated. It has been found that the direct 
current modulation response of many QCLs that have been reported in the literature is over-damped since, 
in contrast to diode lasers, the photon lifetime inside the optical cavity in QCLs is usually the longest 
time constant. The modulation bandwidth is also limited by the inverse photon lifetime. At present, only 
quantum well infrared photodetectors (QWIPs) have a bandwidth wide enough that they could be used to 
study the modulation response of QCLs. However, the current noise provides an alternate way of studying 
the high speed dynamics of QCLs, and as shown in this paper, the modulation bandwidth of QCLs can be 
found by looking at the spectral density of the current noise in the external circuit. 

Although in this paper the emphasis has been on a specific multiple quantum well QCL structure, the 
theoretical methods and techniques presented in this paper can be used to study a variety of QCLs that 
have been reported in the literature. 

Acknowledgments 

This work was supported by DARPA, Rome Laboratories, and ONR. 



A Correlations Among the Langevin Noise Sources 
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B Differential Resistance of a QCL 

The expression in Equations (p4|) and (pfl) for the differential resistance Rd of a QCL can be put in the form, 
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The dimensionless parameters 63, 9 3 , 62, 9' 2 , and 6\ that have been used in the above Equation are as follows, 
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C Elements of Matrix D 

The non-zero elements of the matrix D are, 
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D Important Elements of Matrix D 1 

Above threshold, elements of the matrix D" 1 in the limit u>Ti n <C 1, and {t2 , t{\ — ► oo are as follows: 
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where H(w) is, 
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In the limit WTi n <C 1, and {t 2 , Ti} — > 00 the current modulation response function H{ui) has only three 
poles (Equation (162)). In order to approximate H{lo) with a function with only two poles, as in Equation 



(163), one may neglect the pole which has the largest magnitude. This approximation will be valid provided 
that the pole which is dropped has a magnitude much larger than the other two poles. In this paper we have 
simply dropped the term cubic in u> in the denominator of H(lo) in order to obtain a two-pole approximation 
for H(uj). Numerical simulations show that this is an excellent approximation. In this approximation uir 
and 7 are, 
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E Noise Spectral Densities and Fano Factors 

The spectral densities Ki(uj) and Kp(uj) of noise power for the current noise and the intensity noise, respec- 
tively, can be computed from the equations, 
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Equations (|166| ) and ( |167| ) can be used with Equations (|98|) and (|96|) to compute the noise spectral densities. 
Since all the Langevin noise sources are delta correlated in time domain, they will also be delta correlated 
in frequency domain, and therefore the fluctuations SI and SP ou t in the current and the output power, 
respectively, will be also be delta correlated in time and frequency domains. 

The Fano Factors Fi(uj) and Fp(lo) for the current noise and the intensity noise, respectively, are defined 
as the ratios of the actual noise spectral densities to the noise spectral densities of shot noise, and are given 
by the relations, 
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The Relative Intensity Noise (RIN) is defined as, 
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F Current Noise and the Equivalent Circuit Models 

Expression for 8P{uj) in Equation (|97|) can be found from Equation j95|) and Equations (|S9|)-(|92|), 
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Figure 24: Active region of a semiconductor quantum well diode laser. 



G Noise Model for Semiconductor Quantum Well Diode Lasers 



A simple model for the current and intensity noise in quantum well interband semiconductor diode lasers 
is presented pi, 



24 



The carriers are 



i0{ . The active region of a quantum well diode laser is shown in Fig. 
injected from the leads into the separate confinement hcterostructure (SCH) region cither by tunneling or 
by thermionic emission over the hetero-barrier. The rate equations for the fluctuations SN C and SN W in the 
carrier densities (cm -3 ) in the SCH region and the quantum wells, respectively, and the fluctuations SS P in 
the photon density (cm -3 ) are || |3p|] , 
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It is assumed that the carrier density N c in the SCH region also includes the carriers inside the quantum 
well barriers and also those in the quantum wells which have energy high enough to not be confined within 
the quantum wells (Fig. pi]) . Only those carriers which are confined within the quantum wells are included 
in the carrier density N w . V c and V w are the volumes of the SCH region and the quantum wells, respectively. 
Vp is the volume of the optical mode. r c and r e are the capture and emission times for electrons going into 
and coming out of the quantum wells, respectively. 77 is the lifetime associated with carrier leakage and 
recombination in the SCH region. r nr is the non-radiative recombination time in the quantum wells. r sf 
is the differential lifetime associated with stimulated and spontaneous emission into the lasing mode. r p is 
the photon lifetime inside the laser cavity. SVd is the fluctuation in the voltage across the active region. 
The conductance G relates the increase in the injection current into the SCH region from the leads with 
the increase in the voltage across the active region at a fixed carrier density. t g relates the decrease in the 
current injection rate to the increase in the carrier density in the SCH region. F in is the Langevin noise 
source associated with carrier injection into the SCH region. Fi, F c , and F e model the noise in carrier 



leakage, carrier capture and carrier emission events. F nr describes the noise in non-radiative recombination 
in the quantum wells including spontaneous emission into the non-lasing modes. Frn and F R s model the 
noise associated with photon emission into the lasing mode. Fl and F a model the noise in photon loss from 
the cavity. All the non-zero correlations of the Langevin noise sources can be obtained from the methods 
described in Ref. M, 
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Fi n has the approximate correlation |3l| , 
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The inclusion of the rate equation for fluctuations in the carrier density in the SCH region is necessary to 
accurately model the current noise. Carrier leakage in the SCH region results in a less than unity efficiency 
rji for current injection into the quantum wells, 
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and above threshold the expression for the output power can be written as, 



Pout = VoVi— {I ~ hh) 

q 



(188) 



(189) 



where r\ Q is the output coupling efficiency 



G.l Modulation Response 

The current modulation response of diode lasers follows from the rate equations, and for frequencies less 
than the inverse of the carrier capture time r c , it can be put in the form [EOJ, 
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The relaxation oscillation frequency ojr and the damping constant 7 are, 
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In diode lasers 7 is small at threshold, and the 3 dB frequency ^3 dB increases with the bias current until 
7/v^ equals lor. As the bias current is increased beyond this point the modulation response becomes over- 
damped, and W3dB starts to decrease. Thus, in diode lasers the value for cj 3 dB | max can be found by equating 
7/1/2 with u)r, and it comes out to be, 
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G.2 Differential Resistance 

The differential resistance of the laser diode below and above threshold can be derived from the rate equations 
by removing all the noise sources, 
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where the emission efficiency r/ e is r nr /(r e + r nr ). Values of the time constants r c , r e , r 9 , ti and r nr are 
typically 10 ps, 40 ps, 50 ps, 60 ps, and 1 ns, respectively [f33| . It follows that 9' and 9 have the values 0.97 
and 0.17, respectively. 

The discontinuity ARd in the differential resistance at threshold becomes, 
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Below threshold the current- voltage characteristics of a laser diode resemble that of an ideal PN junction [p2| , 
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where m is the diode ideality factor with values usually between 1.5 and 2. Therefore 
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The above equation shows that the discontinuity in the differential resistance at threshold is K B T/qI t h times 
a factor which is close to unity. 
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Figure 25: Circuit model for the current fluctutions in semiconductor diode lasers 



G.3 Differential Impedance 

The differential Impedance Z(u) of a diode laser above threshold can be expressed in terms of the modulation 
response H(u>), 
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The differential resistance Rd above threshold, given by Equation (195), equals Z(u> = 0). 



G.4 Current Noise 

As in the case of QCLs, the fluctuations 8I(u>) produced by the current noise source that sits in parallel 
with the laser diode (Fig. 25) can be found by looking at the current noise in the external circuit when the 
voltage fluctuation SVd across the diode is zero (because all external sources and impedances are assumed 
to be shorted). This implies that, 



Below threshold, and for frequencies less than W3dB, 5I(u>) is, 
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Equations (203) and ( pO40 show that above threshold the noise associated with carrier injection into the 
active region is not suppressed (since 6 C 1). Above threshold, the carrier density in the quantum wells 
is strongly damped, and only the carrier density in the SCH region provides negative feedback to suppress 
the noise associated with carrier injection. Below threshold, the carrier injection noise is suppressed (since 
9' w 1). Below threshold, the carrier densities in both the SCH region and the quantum wells provide 
feedback to suppress th e ca rrier i nject ion noise. The spectral density Ki(lu) of the current noise 8I(uj) 
follows from Equations (|203| ) and ( |204[ ) , 



(I < Ith) 



W < ^3 dB 



ql + 2ql t 



Hi 



(1 + f) 



+ 2qn sp rj % (I - I th ) 



(i + ef 



(i > ith) 
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In the limit lo — > oo, the current noise is just the noise associated with carrier injection into the active region 
and has the spectral density, 



Ki{u>) 



ql (1 + 26') (Kith) 

ql (1 + 20) +2q I th (0'-e) (I>Ith) 



(206) 



G.5 Suppression of the Current Noise by Large External Impedance 



The current noise SI ex t(uj) in the external circuit in the presence of an external impedance Z s (lo) and an 
external voltage noise source SV S (ui) is (Fig. p5| ), 



SI ext (u>) 



SV s {lo) 



Z(u) 



(Z{lu) + Z s (u)) (Z(lj) + Z.(u)) 



6I(u) 



(207) 



where Z(u>) is the differential impedance of the active region and Z(u> = 0) = Rd- The external impedance 
Z s (tjj) is the Thevenin equivalent of the external circuit impedance and the impedance associated with 
the laser parasitics (ohmic contact resistance, depletion layer capacitance e.t.c). Assuming that SV s (lu) 
represents only the thermal noise originating in Z s (u>), the spectral density Ki ext (u>) of the current noise in 
the external circuit becomes, 



KvM 



\Z{u) + Z„{u)Y 



Z(u>) 



Z{u) + Z s {lo) 



2K B TReal{Z s (uj)} 

\z(u) + z s {uj)Y 



Z{u) 



Z{uj) + Z s {uj) 



K I\ U ) = „2 1" ry,, — ry I, _\ M 
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When Z s (u>) is much larger than the differential impedance Z(u>) of the active region then the current noise 
in the external circuit is just the thermal noise originating in the impedance Z s (u). When Z s (lo) is much 
smaller than Z{u) then the current noise in the external circuit is the noise originating inside the active 
region. By making the impedance Z s (lo) very large the current noise in the external circuit can be suppressed 
well below the shot noise value. 



G.6 Intensity Noise 

Above threshold, and for frequencies less than lu^^Bj SP out (u)) is, 

5P out ((j) = T) i]i hu Iext ^ + rj hv\ (1 - r/j) [V c F c (uj) - V w F e (uj)] ~ rnV c Fi(u)) 



V w F nr {io) - V w F rn (lu) + 1Tst [V p F RS (iu) - V p F L (cu)} \ + F {lo) 



(209) 

High impedance suppression of the current noise 5I ex t(io) in the external circuit can have a profound 
effect on the laser intensity noise through the first term on the right hand side of the above Equation. If 
5I e xt(ijj) is suppressed then the spectral density Kp(uj) of the intensity noise, for frequencies less than ^3dB, 
is, 



hvP ou 



i , 1 . 1 -Vi 
I - i] + 2t] n sp [ 1 — 



w < uj 3 dB 

+ (Vo hvf 

The low frequency Fano factor of the intensity noise at large bias currents becomes 



mi 1 - Vi)- + m — 
? q 



F P {w) 



((J < W 3 dB , I > Ith) 
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(211) 



In diode lasers both r\ Q and r\i have typical values around 0.85, and therefore high impedance suppression of 
the current noise in the external circuit can result in more than 5 dB suppression of the laser intensity noise 
below the shot noise value. On the other hand if the external impedance Z s (u>) is much smaller than the 
impedance Z(uj) of the active region then it can be shown that at large bias currents the Fano factor of the 
laser intensity noise approaches unity, 



F p {lo) 



(lo < CJ 3dB , I > hh) 



1 



(212) 



In practice it is difficult to make Z s (w) very small. One way of obtaining a small external impedance Z s (u>) 
is by using the circuit B shown in Fig. [?] and shorting the RF port of the bias T. At frequencies of interest 
Z s (lo) would then just be the parasitic impedance associated with the laser device and would be dominated 
by the resistance of the device ohmic contacts. 



Table 1: Device Parameters Used in Numerical Simulations (From Ref. |DJ) 



Parameter 


Value 


Lasing wavelength A 


5 . /m 


Operating temperature 


20 K 


Number of gain stages N 


25 (unless stated otherwise) 


Total confinement factor y^ v , T- 7 


er/(0.019JV) w 0.02 AT 


Cavity width W 


11.7 /im 


Cavity length L 

J o 


3 mm 


Facet reflectivities n , T2 


0.27 


Cavity internal loss on 


11 cm" 1 


Mode effective index n F tf 

W J 


3.29 


Mode group index n g 


3.4 


Differential gain a 


- 4.0 x 10~ a cm 


Tin > ^~out > ^"3 


1.0 ps 


T2, Tl 


00 




2.1 ps 


T31 


3.4 ps 


T21 


0.5 ps 




0.31 /iF/cm 2 


C3 > C2 


0.56 /iF/cm 2 


Cx 


0.81 /iF/cm 2 


Xin > Xout 


~ 1 
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